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Abstract 



The first observation of the decay K + — > 7r + 77 is reported. A total of 31 
events was observed with an estimated background of 5.1±3.3 events in the 
7r + momentum range from 100 MeV/c to 180 MeV/c. The corresponding 
partial branching ratio, B(K + -> 7T+77, 100 MeV/c < P n + < 180 MeV/c), is 
(6.0±1.5{stat}±0.7{sys» x 10" 7 . No K+ — > 7r + 77 decay was observed in the 
7r + momentum region greater than 215 MeV/c. The observed tt + momentum 
spectrum is compared with the predictions of chiral perturbation theory. 

PACS numbers: 13.20.Eb, 12.39.Fe 
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We report the first observation of the rare decay K + — > 77+77. This decay provides 
a stringent test of chiral perturbation theory (ChPT) Jl|] since there is no tree-level 0(p 2 ) 
contribution and the leading contributions start at 0(p 4 ) 0. It includes one undetermined 
coupling constant, c, at 0(p 4 ) for which various models predict different values of order 
of 1 (0(1)) Both the branching ratio and the spectrum shape of K + — > 7r + 77 are 

sensitive to c, and the total branching ratio is predicted to be (0.4 — 1) x for c of 
0(1). Measurements of Kl —>■ 77°77 |7],^] indicate the necessity for next-to-leading order 
(0(p 6 )) "unitarity" corrections flJTOl which can be deduced from an empirical fit of the 
decay amplitude of Kl — > 77+77~77°. Similar corrections for K + — > 7T + 77 from K + — > 
7T + 7r _ 7r + predict a 30-40% higher branching ratio than the 0(p A ) calculations with c ~ 
0(1) [0]- From the measurement of the branching ratio and spectrum shape of K + — > 
7r + 77, it is possible to determine a value of c and also to examine whether the unitarity 
corrections are necessary. Observation of this decay mode in excess of the Standard Model 
rate would suggest the existence of new phenomena such as sequential decays of the form 
K + — > 7i + X°, X° — > 77, where X° is a massive short-lived neutral scalar particle. Our 
previous search for K + — > 7r + 77 was confined to the region of phase space with P n + > 
215 MeV/c. We obtained a 90% C.L. upper limit on B(K + -> 77+77) of 10~ 6 when a 
constant matrix element (phase-space distribution) was assumed, but only 10 -3 -10™ 4 when 
the ChPT spectrum with c ~ 0(1) was assumed. The present work expands the search to 
include the region P n + < 180 MeV/c which is favored by ChPT. 

The E787 setup at the Brookhaven Alternating Gradient Synchrotron has been described 
in detail in Ref. [|l^]. An 800-MeV/c kaon beam was slowed by a BeO degrader and stopped 



in a segmented scintillating fiber target located at the center of the detector. Charged 
decay products were momentum-analyzed in a cylindrical drift chamber placed in a 1-T 
solenoidal magnetic field. Their kinetic energy and range were measured by using the target 
and a 15-layer plastic scintillator (range stack) which surrounded the drift chamber. Signals 
from the range stack were fed to 500 MHz transient digitizers (TD) in order to identify 
the 7r + — ► n + — > e+ decay chain in the stopping layer. Photons were detected in a lead- 
scintillator barrel calorimeter surrounding the range stack. The barrel calorimeter had 48- 
fold azimuthal and 4-fold radial segmentation with a total thickness of 14 radiation lengths. 
It covered a solid angle of about 377 sr. Most of the remaining solid angle in the end regions 
was covered by 12-radiation-length end-cap calorimeters located upstream and downstream 
of the drift chamber. 

To avoid background from the two-body decay K + — > 77+77 (K^) with a 21% branching 
ratio and a monochromatic 77+ momentum of 205 MeV/c (K n 2 peak), we searched for the 
K + — > 7r + 77 decay in the regions above the K n2 peak, P n + > 215 MeV/c (77771 region), 
and below the peak, 100 MeV/c < P w + < 180 MeV/c (77772 region). To select K + —> 77+77 
candidates, multilevel triggers followed a threefold strategy: (1) a K + must stop in the 
target; (2) the decay particle must penetrate into the range stack with the range longer 
(77771) or shorter (77772) than that of a 77 + from K^-, and be identified as a 77 + ; and (3) the 
number of photon clusters in the barrel calorimeter must be exactly two with two-photon 
invariant mass (m 77 ) smaller (7777I) or larger (77772) than the 77° mass. The 7777I trigger 
which satisfied the above strategy was achieved by modifying the logic described in Ref. 



12] to reject events with photons showering in the end-cap calorimeters and events with 



m 77 greater than 130 MeV/c 2 (P w + = 207 MeV/c). For the search in the 77772 region, 
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requirements in addition to the beam and range conditions at the first trigger level were (i) 
a large energy deposit (~ 190 MeV) in the barrel calorimeter and (ii) the large 7-7 opening 
angle (# 77 > 75 deg) so as to reduce K n2 background. At the second trigger level, the on-line 
TD requirement to identify the 7r + — > /i + z/ M decay (T^j was removed to avoid an unnecessary 
acceptance loss. At the third trigger level, m 77 was calculated, and an event was rejected 
if m 77 < 190 MeV/c 2 (P n + = 183 MeV/c). From the data taken in 1991, the total number 
of stopped kaons was 3.1 x 10 10 for the 71771 trigger and 6.1 x 10 10 for the 77772 trigger. 
The difference between the two K + counts was mainly due to a hardware prescale factor of 
2 applied to the 71771 trigger. A total of 7.3 x 10 5 events survived the 7T771 trigger, and 
2.7 x 10 6 events survived the 7T772 trigger. 



For the off-line search in the ^771 region fl41| , we required unambiguous identification 
of a single charged track as a tt + by measuring its range, energy and momentum, and 
by observing the 7r + — > fi + decay. Adjacent hit modules in the barrel calorimeter were 
grouped to identify two isolated photon showers (clusters). The hit position in each module 
along the beam axis (z) was calculated from the end-to-end time and energy differences. 
The azimuthal angle (0) of the hit position was determined up to the segmentation of the 
modules. Then, the location of the photon shower in z and <ft was obtained by an energy- 
weighted average of the hit positions. Events were rejected if the photon-pair opening angle 
and an energy of either the high or low energy photon were consistent with K + — > 7r + 7r°, 
7i° — > 77 kinematics. Cuts were applied on the photon energy ratio and opening angle to 
remove events which could be affected by fluctuations in visible energy or by accidentals. 
Finally, the remaining events were tested by a constrained kinematic fit for consistency 
with K + — > 7r + 77. The constraints included total momentum and energy conservation and 
consistency of the charged track's energy, momentum and range with a pion hypothesis. 
Events that had a \ 2 confidence level, P(x 2 ), greater than 0.1 were accepted. No event was 
seen in the region of E n + > 117 MeV/c 2 , or equivalently P n + > 215 MeV/c. 

The analysis for the 7T772 region followed procedures mostly similar to those of the 
7T771 analysis with the following differences. The major background sources in the 7T772 
region are the K + — ► 7r + 7r°7r° (K n3 ) and K + — ► 7r + 7r°7 (K n2 -y) decays with escaping or 
overlapping photons; these channels do not involve muons in the decay products. This 
allowed removal of the 7r + — > fi + tagging condition in order to increase the acceptance by 
more than a factor of 2. Since there was a large correlation between the measured energy and 
range for a low energy 7r + , the range- momentum relation was removed from the constraints 
of the kinematic fit. Photon vetoes were tightened to reject and K n2 -y events with 
escaping photons. A cut was placed on the angle between the tt + and the higher-energy 
photon (^7r+ 71 ) to remove events with a highly asymmetric photon-pair decay. Events were 
rejected if cos6 l 7r + 7l < —0.8. This cut also reduced K n2 ^ events with the radiative photon 
overlaying one of the two photons from tt° decay. Since the radiative photon from K n2l was 
preferentially aligned along the the ir + direction, an overlapping-photon event usually has 
one 7T° photon near the 7r + direction and the other in the opposite direction. Events were 
accepted if the constrained kinematic fit satisfied P(x 2 ) > 0.1. At this point, the background 
with a photon escape was estimated to be less than one event. 

Remaining background events were primarily due to K n3 and K n2l decays with multiple 
photons fused in the barrel calorimeter to form a single cluster. Background with overlapping 
clusters was greatly reduced by rejecting events with cluster shapes inconsistent with that 
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of a single photon. A typical single photon cluster contained 4 - 5 hit modules in the barrel 
calorimeter. The maximum discrepancy among z measurements and the standard deviation 
of 4> measurements were calculated for each photon cluster, and events were rejected if one 
of the cut variables exceeded the 72% acceptance point. The acceptance was measured using 
a clean single-photon cluster based on K n 2 data with the same number of hit elements and 
a similar energy. A total of 31 K + — ► 7r + 77 decay candidates survived all the cuts (Fig. [[]). 
The residual background was estimated by measuring the rejections of the z and <p cuts. 
Each cut was inverted to enhance the overlapping-cluster background, and the rejection of 
the other cut was evaluated by assuming the two cuts were independent. The independence 
of the cuts and the magnitude of the background level were confirmed by a large statistics 
Monte Carlo simulation. The estimated number of the total background was 5.1±3.3. 

The acceptance was estimated by using Monte Carlo simulation and calibration data 
sets. The net acceptance for the ^771 trigger and analysis was (1.5 ± 0.2) x 10~ 4 which 
included a factor 0.12 for the 7r + momentum being between 215 MeV/c and the end point 
227 MeV/c assuming the phase-space distribution. Under this assumption, the 90% C.L. 
upper limit was B(K + — > 7r + 77, phase-space distribution) < 5.0 x 10~ 7 . 

The momentum dependence of the 71772 acceptance is shown in Fig. |l] together with the 
final 7T77 candidates and the estimated background spectrum. A partial branching ratio 
for each 10-MeV/c bin was calculated from the background-subtracted signal divided by 
the number of the stopped kaons times the acceptance for that bin. By summing these, 
a model-independent branching ratio B(K + — > 7r + 77, 100 MeV/c < P 7i + < 180 MeV/c) 
was obtained to be (6.0 ± 1.5 ± 0.7) x 10~ 7 where the first uncertainty is statistical and 
the second is the estimated systematic uncertainty in the background and the acceptance 
measurements based on a K n2 branching ratio measurement. 

The measured 7r + spectrum was compared with ChPT predictions integrated over 10- 
MeV/c bins. A maximum likelihood fit of c to the spectrum using the absolutely-normalized 
rate was carried out to determine the value of c. Without the unitarity corrections, the 
result was c = 1.6 ± 0.6 with xLin. = 6-3 (n.d.f. = 7). With the corrections, the best fit 
was obtained for c = 1.8 ± 0.6, and the Xmin. improved to 4.6. Thus, the data support the 
inclusion of the unitarity corrections. If we fit the spectrum shape alone by introducing 
one free normalization factor, the fits became equally good, with Xmin./ n -d-f- ~ 0.7. The 
corresponding c was —0.61}^ and 0.71^ for without and with the unitarity corrections, 
respectively. Fig. |2| shows the measured momentum spectrum together with the spectrum 
shape given by the best fit (i.e. c = 1.8 with the unitarity corrections). By assuming this 
spectrum shape, a total branching ratio was estimated to be (1.1 ± 0.3 ± 0.1) x 10 -6 with 
the net acceptance in the ^772 region being 3.8 x 10~ 4 . The total branching ratio estimated 
with a spectrum shape given by the other fits was within 11 % of the above value. 

One of the consequences of the unitarity corrections is a non-zero amplitude at the end 
point P n + = 227 MeV/c (m 77 = MeV/c 2 ). However, the predicted decay rate at the 
end point is 8 times smaller than our 90% C.L. upper limit in this region, B(K + — > 7r + 77, 
P^ > 215 MeV/c ) = 6.1 x 10~ 8 , as shown in Fig. £| Therefore, the ChPT prediction with 
the unitarity corrections is consistent with our 71771 result. 

The 7T771 result also sets a 90 % confidence upper limit on B(K + — > 7r + X°,X° — > 77), 
where X° is any short-lived neutral particle with a mass smaller than 100 MeV/c 2 decaying 
into two photons. Fig. ^ shows the limit as a function of M x o for different lifetimes. 
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FIG. 1. 7r + momentum distribution for 31 K + — > 7r + 77 candidates (solid) and for estimated 
background events (dashed). The acceptance is given by the dotted line. 
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FIG. 2. Measured ir + momentum distribution for K + — > 7r + 77 and the best fit to the data 
(solid). The dashed line shows a phase-space distribution normalized to a 90% C.L. upper limit 
obtained by the ^771 analysis in the region of 215 MeV/c < P w + < 227 MeV/c indicated by the 
arrow. 
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FIG. 3. The 90% C.L. upper limits for the branching ratio of K + -> 7r+X°,X -> 77 for 
different X° lifetimes (r^o ) as a function of mass (m x o ) • 
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